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Abstract 

Excitonic spectra of armchair graphene nanoribbons (AGNRs) obtained from a full many-body 
exact diagonalization of the Hubbard model are reported for both longitudinally and transversely 
polarized photons, thus providing a complete survey of low-energy may-body optical excitations in 
these systems. The resulting one-photon allowed eigenstates turn out to be well separated in energy 
from each other but both couple to the same set of two-photon allowed states. The magnitude of 
the calculated optical oscillator strengths for perpendicular polarization suggest that these optical 
features can be indeed observed in polarized absorption measurements. 
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I. INTRODUCTION 



The recent progress in the synthesis and isolation of single graphene-layers [l-3| has con- 
siderably increased the interest in graphene nanoribbons (GNRs), infinitely long stripes of 
carbon atoms arranged in a honeycomb lattice. The possibility of tuning the width of these 
quasi-one- dimensional (ID) nanostructures alongside their edge geometry provides a strat- 
egy to overcome the absence of an electronic gap in graphene, which is one of the main limits 
preventing the application of this material in electronic devices ^, |5j. Recently, high quality 
stripes less than 10 nm wide, which can be exploited for making field-effect transistors, have 

riQ 

been indeed fabricated from bottom-up techniques [6H9|, thereby opening the path for a 
range of experimental investigations on the optical and transport properties of these sys- 
tems. This will be a crucial benchmark for testing the theoretical predictions provided until 
now for GNRs, which have all been recently surveyed in several excellent reviews 13Nl6|. 



Most of these studies rely on independent-particle approximation with proper boundary 
conditions at the GNR edges, such as tight-binding (TB) methods for tt electrons 10|, the 
k ■ p two-dimensional Weyl-Dirac equation for free massless particles whose Fermi velocity 
(lO^m/s) plays the role of an effective speed of light [ll| and ab initio density functional 
calculations |l2|. However, since electron-electron interactions are expected to be enhanced 
in such low dimensional systems, there have been also several theoretical investigations con- 
cerning the inclusion of many-bod y e ffects beyond single-particle theory. These rely on ab 
initio GW-Bethe Salpeter methods|26l. 2?! , the solution of Pariser-Parr-Pople effective model 



Hamiltonians for vr electrons incorporating longer-range Coulomb interactions |28 



29| and 



Hubbard model based approaches, either within mean-field approximation p/7[ or beyond it 
with configuration interaction (CI) carried out at different levels 30, |4l|. The suitability 



of mean-field approximation for the Hubbard model when applied to graphene-based sys- 
tems has been matter of recent debate in literature 



31 



34j . in view of the low- intermediate 



value for the Hubbard correlation coupling strength {1 < U/t < 2.2) which seems plausible 
for these systems. The geometrical classification of GNRs is dictated by their edge shape, 
which can be zigzag (Z), armchair (A) or chiral (C). Moreover, the structure of armchair 
nanoribbons (AGNRs) can be derived from that of zigzag single-walled carbon nanotubes 
(ZSWCNTs) by unzipping the graphene cylinder along the nanotube axis, while conversely 
zigzag GNRs (ZGNRs) can be obtained from armchair SWNTs. This concept has been ac- 



tually exploited for producing controlled-sized GNRs through chemical jsl, Q and laser 36 1 
unzipping of multi-walled CNTs and it has also been recently extended to the fabrication of 
graphene quantum dots starting from fullerene cages \3T\ \. 

The absence of zero-energy localized states simplifies considerably the investigation of the 
electronic properties of AGNRs in comparison with ZGNRs. Besides that, the electronic 



band picture provided for AGNRs closely resembles that o: 



ZSWCNTs, since the van Hove 



singularities occur at the centre of the Brillouin zone |17H19|. On the same footing, the 
optical selection rule for AGNRs with light polarization parallel to the GNR edge is clearly 
reminiscent of that found for conserved-quantum number interband transitions in ZSWC- 
NTs, as discussed in several recent works reporting on optical selection rules and analytical 
expressions for the TB electron-light interaction matrix elements in both armchair and zigzag 



GNRs with either longitudinally or transversely polarized photons 



22|-|25|. 



Both TB and Weyl-Dirac equation predict that armchair GNRs with pristine edges may 



owever, 



20| and 



be either semiconducting or metallic according to their width (oscillating gap). H< 
ab initio and mean-field Hubbard model results, recalled in the reviews by Cresti 
Rozhkov 21[, point out an always semiconducting behaviour for AGNRs, since the metallic 
state is unstable against bond deformations at the edges, electron-electron interactions and 
longer-range hoppings. 

This work expands our previous full-many body exact diagonalization (ED) results for the 
Hubbard model a ppli ed to investigate optical excitations in AGNRs with longitudinally 
polarized photons 4l| by considering also the polarization component perpendicular to the 
GNR edge, thus providing a complete overview of excitonic effects in these systems. Two-leg 
ladder models mimicking AGNRs of different widths and belonging to three distinct families 
are used to sample the set of k-points lying at the centre of the Brillouin zone which give 
the van Hove singularities and the maximum absolute values of the interband optical matrix 
elements. The appropriate velocity operator for optical transition is mapped onto these 
quantum lattice models to give the chosen polarization component in the calculation of one- 
and two-photon absorption optical matrix elements. The obtained trends for the allowed 
transitions and related optical oscillator strength as a function of the correlation coupling 
strength U/t provide an insight into the family dependence of the optical anisotropy in AG- 
NRs and suggest that the spectral features related to transversely polarized photons can be 
observed in these systems in polarized absorption measurements. 
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II. METHOD 



In FigJTJwe show the geometrical structure of an armchair graphene nanoribbon unit celL 
Its width W is equal to the number N of dimer lines containing A-type and B-type carbon 
atoms, hence a total number of 2N sites. Periodic boundary conditions along the direction 
parallel to the edges are represented through dashed lines. Following the traditional nomen- 



clature for GNRs 



12| . we consider in this work AGNRs with 4 < < 7. Incidentally, we 



note that such small-width structures are not merely hypothetical, since the 7-AGNR struc- 
ture has recently been obtained from aromatic precursors and investigated both by scanning 
tunneling microscopy and Raman spectroscopy [9]. 

Nth dimer line — > 
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FIG. 1. (color online) Armchair nanoribbon of width W with N dimer lines. ^(i?)-type atoms are 
shown in red (blue), respectively. Equivalent atoms on the honeycomb lattice are mapped by the 
ai and 02 basis vectors. 



Since A''-AGNRs are topologically equivalent to to brick type lattices (i.e. with periodic lad- 



ders) [l2|,|20| and periodic boundary conditions along the direction parallel to the nanoribbon 



edges are considered {k\\ =0 states are sampled in the GNR Brillouin zone), one can fold 
the brick type lattice into a two-leg open ladder having rungs and a unique value t for 
all nearest neighbour hopping parameters in the Hamiltonian (see Fig. 3 in Ref. [3]). We 
consider the simple Hubbard Hamiltonian form for vr electrons 



(1) 



where i and j are site indices, are all pairs of first nearest neighbor sites, cj^ and c^^ 
are electron creation and annihilation operators, rii^a- = c\ ^c^ ^ is the number of electrons on 
site i with spin cr, and U are the nearest-neighbour hopping parameter and the on-site 
Coulomb repulsion parameter between two electrons with opposite spins, respectively. We 
recall that within this Hamiltonian the on-site Coulomb interaction U must be considered 
an effective parameter, whose value for half-filled systems with non-small values of the 
interactions can be taken equivalent to V^^^ = U — Vi, where is the first nearest neigh- 



bour Coulomb interaction in the extended Hubbard models [29|, l39|]. Thus the (effective) U 
parameter takes implicitly into account also longer range Coulomb interactions in the limit 
of static screening, as discussed by Wehling et al. [38]. 

The one-photon optical spectral function of the AGNR is calculated according to the 
Lehmann representation 

I{E) = 5^ |(^„|t;^^|^G5)|'5 {E + Egs - E^) (2) 

m 

where Egs is the ground-state (GS) energy of the system and Em the energy of any other 
eigenstate {ipm) obtained from exact diagonalization of the Hubbard Hamiltonian and v"^^ 
is the velocity operator for light polarization either along (||) or perpendicular (±) to the 
GNR edges, whose general form can be expressed concisely as 

''^"--xE (<=U»-4.U (3) 



The extended expressions of f^^y are reported for completeness in the Appendix for the 
4-AGNR sample system and can be straightforwardly derived for the remaining structures. 
The operator t;^^ obeys the selection rule for the azimuthal band quantum number /i, namely 
A/i = for interband transitions in AGNRs with light polarization along the nanoribbon 



edges 



22 



2^. As it occurs with zigzag SWCNTs 



39 



401 1 , this implies that vertical transitions 



occur between van Hove singularities belonging to the same AGNR band quantum number 



in the TB picture, as noted in our previous work on AGNRs 



411 ]. There we made the choice 



of considering only the case of polarization along the ribbon edges as in most literature 
on excitonic effects in AGNRs, since in quasi-lD materials strong depolarization effects 
are believed to quench the optical absorption for perpendicular polarization j^^]. However, 
because of the structural anisotropy of AGNRs, we expect that the optical response of these 
systems should exhibit anisotropic features which could in principle be detected in polarized 



absorption measurements, also observed by Gundra and Shukla 28|]. In order to test this 
hypothesis, we determined the optically allowed transitions also for the operator, which 
takes into account vertical transitions with non-conserved band quantum number, that is 
with Afi = ±1,±3,... 

The non-zero matrix elements of v"^^ between each generic eigenstate |m) and the one- 
photon allowed eigenstates previously determined through Eqj2] allow us to recognize the 
two-photon allowed states in the obtained eigenset. 

We consider a half-filled system with n electrons distributed over n = 2N sites and total 
spin quantum number 3^ = 0. Thus the size of the basis set and the dimension of the 
matrix to be diagonalized is D = [n\ / {n-t-ln^l)]'^ , where and ni are the numbers of spin-up 
and spin-down electrons, respectively, with n-^ = = N. Calculations were performed for 
several values of the U/t correlation coupling strength. An additional set of calculations 
was performed for non-zero Sz, namely Sz = {n-^ — ni)/2, in order to identify among the 
obtained Sz = eigenstates the corresponding spin multiplicity. In this way, only those 
Sz = eigenstates which also appear in the Sz = I set but not in the Sz = 2 set are 
classified as triplet excitations, whereas the remaining Sz = eigenstates, which do not 
belong to the Sz = 1 set at all, are classified as singlets. We verified that in the low-energy 
spectrum, it is sufficient to perform ED up to the Sz = 2 quantum number, since states 
belonging to the Sz = 3 set do occur well above the energies of the states of interest for the 
low-energy optical properties. In order to perform ED for the Hubbard model, we adopted 
an iterative diagonalization scheme based on the Lanczos algorithm, as implemented in 
the ALPS libraries I42I, and additional matrix- free strategy combined with shared- memory 
parallelization on multicore nodes 43| in order to speed-up the diagonalization of the larger 
systems with 12 and 14 sites mimicking the 6- and 7-AGNRs, respectively. 



III. RESULTS AND DISCUSSION 

In Figj2] a general overview of the one- and two-photon allowed transitions is given for 
either longitudinally and transversely polarized photons in AGNRs with = 4 7 in 
the correlation regime < U/t < 4. One-photon allowed transitions obtained from the 
perpendicularly-polarized velocity operator vj^ are denoted by the (+) symbol, whereas the 
ESll and ES22 obtained from the v?P operator are denoted by the (*) and (x) symbols. 
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respectively. We find that both components of the velocity operator v"^^ couple the obtained 
one-photon allowed states to the same set of two-photon allowed transitions. For the case 
U/t = the TB results were verified for the interband optical transition energies and the 
metallicity condition, N = 3p + 2, with p an integer, for 5-AGNR. For the remaining ribbon 
families with N = 3p (6-AGNR) and N = 3p + 1 (4- and 7-AGNR) we verified the semi- 
conducting behaviour. However, we recall that when electronic correlations are taken into 
account U/t > the metallic behaviour of the 5-AGNR is no longer observed because of the 
lifting of the K-point degeneracy and all the three AGNRs families display semiconducting 
behaviour. 

As a feature common to all the three AGNR families, the bright states given by vj^ are 
always found between the ESll and ES22 transitions obtained from f Using the values 



tyr = 2.6 -ir 2.8 eV for the TB hopping parameter in graphene-based materials |31l . |38| . the 
magnitudes of the one-photon transition energies are of the same order of those previously 
computed by GW Bethe-Salpeter methods 0, H] or by the PPP Hamiltonian 28, 29] for 
the corresponding families. 




FIG. 2. Transition energies vs U/t for the considered A'^-AGNRs with N = 4,5,6,7. Symbol 
legends: (*) one-photon first bright state Ell obtained from tiy^; (+) one-photon bright state(s) 
obtained from v'j^; (x) one-photon second bright state E22 from v"^^^; two-photon states are 
denoted by triangles. 
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Interestingly, we note that v"^ activates only one transition in the = 3p, 3j9 + 1 families, 
whereas two one-photon allowed transitions are found for perpendicular light polarization 
for the N = 3p + 2 family (5-AGNR). This finding again confirms our previous result that 
the opto-electronic properties of AGNRs are family- dependent, as also observed by Prezzi 
et al. , although the behaviour of the A^ = 3j9 + 2 family is significantly different from 
those of the N = 3p,3p + 1 systems |41| . 

Since graphene nanoribbons are in the intermediate-low correlation regime 1 < U/t < 2 33| . 
at least one two-photon allowed transition is always found between ESll and ES22 for all 
the three AGNR families, while the one-photon transition given by vY' can be found either 
below {U/t ^ 1) or above {U/t ~ 2) the two-photon allowed state. 

In Fig. [3] the magnitude trends of the optical oscillator strengths for the considered one- 
photon allowed transitions are reported as a function of the correlation strength U/t. 




FIG. 3. Oscillator strengths vs U/t for the considered TV- AGNRs with A^ = 4, 5, 6, 7. Symbol 
legends: (*) one-photon first bright state Ell obtained from v"^^^; (+) one-photon bright state(s) 
obtained from v'^; (x) one-photon second bright state E22 from 'Wy^- 

Again for the A^ = 3p + 2 family, we note a strikingly different behaviour than in the 
other systems, since the t^-normalized optical oscillator strengths of the vY' transitions are 
quite low (about 4 times lower than the ESll or ES22 oscillator strengths), whereas in the 
A^ = 3p, 3p -|- 1 families the reported values for the corresponding transition are comparable 
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or even higher (6-AGNR) than the ES22 transition. Since the transitions are well sep- 
arated in energy from those obtained with longitudinally polarized photons, these results 
predict that anisotropic optical features are likely to be observed more distinctly in AGNRs 
belonging to the = 3p, 3p + 1 families. 
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FIG. 4. Exciton level scheme of 6-AGNR (A = 3p) obtained from ED with Ujt = 1.5. Level 
legends: singlets (black), triplets (blue), quintuplets (red). ESll, ES22|| are the one-photon allowed 
transitions obtained from ^ ESllj_ from w^^, TPAl-4 are the two-photon allowed states. 

In Fig. Hlwe report the low-energy excited-state spectrum of 6-AGNR with an explicit char- 
acterization of the spin multiplicity of the optically inactive (dark) exciton states alongside 
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the allowed one- and two-photon transitions for both polarization directions. The 6-AGNR 
was chosen since according to Fig. |3]this system displays high values of the ESll oscillator 
strengths for both polarization components in the considered correlation regime. Anyway, a 
very similar fine structure of the exciton levels can also be obtained for 7-AGNR. We observe 
the presence of deep-lying triplet excited states above the ground state, as reported in SCI 
calculations by Dutta et al.[30|. Moreover a dark singlet is found 81meV above ESlly for 
t = 2.6 eV, whereas two triplets are found 82meV above and 31 meV below ESll^, respec- 
tively. As discussed in the case of SWCNTs |40|, these states can act as a population sink 
for the nearby optically-allowed exciton through a bottleneck mechanism. In particular, the 
spin-orbit coupling due to impurities or defects is expected to enhance the quenching of the 
radiative transition by trapping the radiation in triplet dark states, so that the luminescence 
quantum yield from the bright state would be decreased. 



IV. CONCLUSIONS 



In summary, we have investigated the excitonic structure of pristine small-width AGNRs 
by ED of the Hubbard model for several two-leg ladder models mimicking ribbons with 
4 < < 7 dimer lines and for both components of photon polarizations, either along and 
perpendicular to the ribbon edges. By this technique we are able to investigate the effect of 
electronic correlations involved in direct interband optical transitions by a full-many body 
approach for several values of the Hubbard correlation parameter U/t. The obtained results 
allow us to track down the nanoribbon family dependence of the anisotropic features in 
the optical properties of these systems, which could be detected in polarized absorption 
measurements. 
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Appendix A: Velocity operators for interband optical transitions in armchair GNRs 



In the following we report the explicit second-quantization expressions for the components 
of the velocity operator f^^n j_ along and perpendicular to the edges of the 4-AGNR. The 
site indexing of the creation/annihilation operators refers to the atomic labels reported in 
the related Fig. [51 

cr 

+ 1 {claC2,a + 4,.C4,. + c{.C^,. - h. C.) + (Al) 

= E [4,.C2,. + 4,.C3,. + 4,.C4,. + 4,.C5,.+ (A2) 

+ 4,aC6,a + 4,aC7,a " h- C. 
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